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Abstract 

The polarization properties of RX J1856. 5-3754 are investigated, 
based on the plasma emission model presented in previous works, 
is shown that if the emission of this source is generated by the 
/nchrotron mechanism, the X-ray and the optical emission are 
inearly polarized and the corresponding degrees of linear polariza- 
jon are calculated. The measurement of the polarization patterns 
If this source appears to be most reliable way, to make clear its real 
smission nature. 



Introduction 



(N : 

C/1RX Jl 856.5-3754 (or RXJ1856, hereafter) is one o f the brightest 
ae'arb y isolated neutron stars discovered by ROSAT ( Walter et al. I 
\^9Q) as an X-ray source and a large number of papers have been 
•devoted to explaining its nature. Its soft X-ray spectrum can be 
p^ell fit by the Planckian sp ectrum with a temperature kT™ b ^ 63 + 
"^eV dBurwitz et al. 1 120031) . Thus it is supposed that the emission 
Cfif this source is of thermal nature. However, the existing models 
1 Based on this original assumption face problems. 
t— | The lack of any significant spectral lines in the X-ray spec- 
>^tru m of RXJ1856 argues agains t heavy element atmosphere mod- 



'^ls dBurwitz et al. , 2001 , 120031) . whereas single temperature hy- 



ro gen atmosphere fits over-predict the opt i cal flux by a la r ge fac - 
r dPavlov et all 1 19961 IPons et al. 1 120021: iBurwitz et ail [2003). 



^ As none of the classic models of neutron star atmosphere are able 

f- tp explain the observed X -ray spectrum, which is well fitted by a 



Q^imple blackbody model dBurwitz et al. L 120031) . it has been pro 



Cpose d that the star has no a t mosphere but a condensed mat ter sur 
2001 



This 



*^ace dBurwitz et al. L I200U iTurolla. Zane & Drake 1 12004 
• ^urface might emit a virtually featureless b lackbody-like spec 
)/Hfum (as originally proposed bv lPavlov Since the broad 



^and spectra of this source can not be fitted by a single Planck- 
ian spectrum it is often described by two-temperature black- 
body mode ls ( Pons e t al. I l2002t IPavlov. Zavlin & Sanwal L 12002c 
Burwitz et al. L |2003| ). However, condensation o f surface matter 



requi r es very spe cific conditions to be fulfilled dLai & Salpeter 



1997; Lai , 2001). But even if these conditions are satisfied, the 



formation of a non-uniform distribution of the surface temperature 
(two-temperature blackbody models) still remains unclear. The 
most adequate fits of the spectra give models which assume that 
the star has a thi n hydrogen atmosphere superposed on a condensed 



matter surface (Mot ch et all 120031: IHo et al. LI2007I) . However, the 



origin of such thin hydrogen layers fitting the data is a problematic 



Alternatively, the observational properties of RXJ1856 can 
be explained in the framework of the plasma emission model 



presented in previous pa p ers dChkheidze & Machabeli L 12007 



Chkheidze & Lomiashvili L 120081) . which are based on well- 
developed theory of pulsars. This model suggests that the emis- 
sion from this object is generated by synchrotron radiation, cre- 
ated as the result of acquirement of pitch angles by relativistic 
electrons during the quasi-linear stage of the cyclotron instabil- 
ity. The model gives successful fits for broadband spectra, with- 
out facing the problems typical for the thermal radiation mod- 
els. Considering the case of a nearly aligned rotator, it was 
predict ed that the source should have pulsated with a period of 
~ 1 s (Ch kheidze & Machabeli 1 120071) . However, the posterior 



XMM-Newton observation of RXJ1856 disc overed that its X-ray 



emiss ion pulsates with a period of 7.055 s dTiengo & Mereghetti 



120071) . The latter fac t has been explained in the framew ork of drift 



wave driven model (Chk heidze & Lomiashvili L 120081) . In partic 



ular, the real spin period of the pulsar might differ from the ob- 
servable one, as a consequence of the existence of very low fre- 
quency drift waves in the region of generation of the pulsar emis- 
sion. These waves are not directly observable but only result in 
a periodical change of curvature of the magnetic field lines and, 
hence, a periodical change of the emission direction with a period 
of the drift waves assum ed to be equal to the observable period 



dLomiashvili et al. 



20061) . 

We are not about to reject the existing thermal radiation 
models for RXJ1856 and suppose that the most reliable ar- 
gument for revealing the real emission nature of this source 
will be its study with polarization instruments. Therefore in 
present paper the polarization properties of RXJ1856 are in- 
vestigated, in the framework of the plasma emission model 
developed in previous works dChkheidze & Machabeli L 12007 : 
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Chkheidze & Lomiashvili 1120081) . 

In this paper, we give a description of the emission model in 
(Sect. 2), investigate the emission polarization in (Sect. 3) and 
discuss our results in (Sect. 4). 



2 Emission mechanism 

It is well known that the distribution function of relativistic parti- 
cles is one dimensional at the pulsar surface, because any trans- 
verse momenta {p ± ) of relativistic electrons are lost in a very short 
time (< 10" 20 s) via synchrotron emission in very strong B ~ 10 n G 
magnetic fields. For typical pulsars, the plasma consists of the fol- 
lowing components: the bulk of plasma with an average Lorentz- 
factor y p ^ 10; the tail of the distribution function with y, ^ 10 4 
and th e primary beam with yt, — 10 6 (see Fig.l from Arons 



1981)). Generation of waves is possible during the further motion 



1 



of the relativistic particles along the dipolar mag netic field lines 



if the condition of cyclotron resonance is fulfilled dKazbegi et al. 
1991b : 



0) - k^Vy, - k x U x + COs/jr = 0, 



(1) 



where k^+k]_ = k , k]_ = k x +k r , is the particle velocity along the 
magnetic field, y r is the Lorentz-factor for the resonant particles, 
u x = cVyjrlpwB is the drift velocity of the particles due to curva- 
ture of the field lines, p is the radius of curvature of the field lines 
and o>b = eB/mc is the cyclotron frequency. Here a cylindrical 
coordinate system is chosen, with the x-axis directed perpendic- 
ular to the plane of field line, when r and <p are the radial and az- 
imuthal coordinates. The reason for the generation of such waves is 
an anisotropy of the distribution function. During the quasi-linear 
stage of the instability, a diffusion of particles arises not only along 
but also across the magnetic field lines. The resonant electrons, 
therefore, acquire transverse momenta and, as a result, start to ra- 
diate t hrough the synchrotron mech anism. 

In Chkheidze & Machabeli (2007), it has been assumed that the 
emission of RXJ1856 is generated by the synchrotron mechanism, 
which switches on during the cyclotron instability developed at dis- 
tances r ~ 10 9 cm. The observed X-ray spectrum is a result of radi- 
ation of the primary-beam electrons. Particularly, the synchrotron 
spectrum of the beam electrons with the energy distribution: 



f\\b 



n b 



exp 



{jr ~ Jb) 



(2) 



well matches the measured X-ray spectrum in the energy interval 
(0.26 - 0.9)keV (here y T 10 - is the half-width of the distribution 
function, «/, = B/Pce is the density o f primary beam electrons , 



1969) 



equal to the Goldreich-Julian density (Gol dreich & Ju lian 
and P is the pulsar spin period. The wave excitation leads to a 
redistribution of the resonant particles via quasi-linear diffusion. 
Consequently, by achieving the stationary mode the distribution 
function of the beam electrons takes the following form: 



f\\b' 



r; 4 . 



(3) 



At the same time the radiation density appears to be sufficiently 
high and prevents the domination of the self-absorption processes. 
Synchrotron self-absorption redistributes the emission spectrum 
and, in the energy domain of relatively low frequencies it expe- 
riences drop. Hence, for the energy interval (0.15 - 0.26)keV the 
X-ray emission spectrum h as the form v 5 ^ 2 , which fit s well the ob- 
served one (see Fig.4 from lChkheidze & Machabeli I (|2007)). The 
measured optical spectrum is the result of the synchrotron emission 
of the tail electrons, redistributed due to cyclotron instability: 



f\\t 



(4) 



The resulting theoretical spectrum matches the measured one 
closely. 



3 Emission polarization 

Now let us consider the polarization properties of RXJ1856. For 
this reason, we have to find the Stokes parameters, which, in our 




Figure 1 : k is the emission direction, B - the magnetic field line, v 
and v' - are the electrons' velocities. 



case, are defined as follows dGinzbureUl981l) : 

/(v) = j^f B [yj +* 2 y 2 )y/(y,<A)sinax 

x [(1 + xV)Kl^g v ) +X 1 y 2 K 2 ll ,{g v )\ d X dy, 



(5) 



G(v) 



2tt 



( l + x 2 y 2 )y 2 .f(y> <A) sm a cos ^-x x 



x [(l +xV)K 2 2/ i(gv) -xV^ig,)} d X dy, (6) 



+ * 2 y 2 )y 2 /(y, $) sin a sin 2 X x 



x [( 1 + X 2 7 2 )K 2 2/3 (8 V ) - xYK^igv)] d X dy, (7) 



V(v) 



where 



= v / B [v) ( 1+ /y 2 ) 3/2 yV(%<A)sinffx 

xK l/3 (g v )K 2 /3 (g v )dxdy, 



4 ■ 10 3 By 2 , 



2v c 



(1 + xV) v \ 



(8) 



(9) 



(10) 



Here a is the angle between B and k, iff - is the pitch angle, X is 
the angle between the wave vector and the electron velocity vector 
(see Fig. Q]). Also, y is the Lorentz-factor of emitting electrons, 
f\\(y) and/ ± (i/0 are the distribution functions of the emitting parti- 
cles by parallel and perpendicular momenta, respectively, K2/i{g v ) 
and K\/3(g v ) are the Macdonald functions and angle x defines the 
direction of maximum intensit y of the polarized component on the 
plane of sky. According to IChkheidze & Machabeli ( 2007 ) the 
following condition is fulfilled for the resonant particles: t//oy r » 1 
(where i/'o is the mean value of the pitch angle). In this case, the 
emitting particles' distribution function by thei r pitch angles has 
the following form ( Malov & Machabeli! 120021) : 



/±(<A) x e 



4e b 



3r:h 



rl 



(ii) 



2 




Figure 2: Point B corresponds to the magnetic field line and point 
A - to the observer's line of sight. 



Synchrotron emission of a single electron is strongly beamed 
along the direction of motion into an angle of order l/y. Hence, in 
the given direction, the observer will detect radiation of the elec- 
trons with velocities filling up the cone with l/y in angular size 
and with the major axis coincided with the line of sight of an ob- 
server (see Fig. [TJ. Therefore, at the given moment of time, the 
observer receives emission of the electrons having the pitch-angles 
from the interval iff = a+x C [a - l/y, a+ l/y] (it should be men- 
tioned that the angle x is considered to be positive \f ip > a and to 
be negative, otherwise). The Stokes parameter V is an odd func- 
tion of variable^, therefore, its integration over symmetric interval 
yields zero. Let us consider the top view of Fig. Q] (see Fig. |2j. 
The larger circle represents a cross section of the cone. The lateral 
surface of this cone is described by the electron velocity vector 
moving on a spiral along the magnetic field line. The circle with 
smaller radius corresponds to the cone combined by the velocity 
vectors giving the significant radiation in the observer's direction. 
Its angular size estimates as l/y. The area of the smaller circle is 
striped in two different ways: area a corresponds to the electrons 
which have positive values of angle x, while area b corresponds 
to those having negative values of angle x- As these areas are not 
equal to each other, one would expect the nonzero value of V. But 
bearing in mind that the bigger circle is large enough in compari- 
son with the smaller one (the solid angles of corresponding cones 
are of the order of 10~ 3 and 10~ 6 for the beam electrons and, 10 _1 
and 10~ 4 for the tail electrons, respectively), one can easily assume 
that a ~ b. But only symmetry of integration bounds is not enough 
condition to have a zero circular polarization. Also the distribution 
function by pitch angles, containing the variable x, should change 
very slowly within the small interval ~ l/y. The estimations show 
that the latter condition is well fulfilled. Thus, we conclude that 
the emission of this source is not circularly polarized. 

Now we can calculate the d egree for the linear polarization, 



which writes as dGinzburgU l981): 



n 



Ve 2 + u 2 



(12) 



The degree of linear polarization, corresponding to the energy in- 
terval (0.26 -0.9)ke V of the observed X-ray spectrum can be calcu- 
lated using Eqs. (2), (5), (6), (7) and (11). We get II; * (84 - 79)% 
for the energies from the given interval. As it was mentioned 
above, the measured X-ray spectrum at relatively low frequencies 
(0.15 - 0.26)keV is produced due to synchrotron self-absorption 
process, of the beam electrons with the distribution f« b oc y ; 7 4 . The 



degree of linear polarization in this case is defined as follows: 

' 6c/+13 



(13) 



where d — 4. The calculations give 11^ « 8%. So, in the energy 
interval (0.15 - 0.26keV), the degree of linear polarization of the 
X-ray emission undergoes the strong drop. 

The observed optical spectrum is the result of the synchrotron 
emission of tail electrons, redistributed due to quasi-linear diffu- 
sion: 

h - r 2 - (14) 

The degree of polarization of the optical emission yields Tlj w 
86%, in this case. 



4 Discussion 



In previous works (Chkheidze & Mac habeli 



2007 



Chkheidze & Lomiashvili , 2008b the plasma emission model 
of RXJ1856 is presented. It is supposed that the emission of 
this source is generated by the synchrotron mechanism. The 
distribution function of relativistic particles is one dimensional 
at the pulsar surface, but plasma with anisotropic distribution 
function is unstable that can lead to a wave excitation. The 
main mechanism of wave generation in plasmas of the pulsar 
magnetosphere is the cyclotron instability. During the quasi-linear 
stage of the instability, a diffusion of particles arises as along, also 
across the magnetic field lines. Therefore, plasma particles acquire 
transverse momenta and, as a result, the synchrotron mechanism 
is switched on. The measured X-ray and optical spectra are the 
results of the synchrotron emission of primary-beam and tail 
electrons, respectively. The predictable characteristic frequencies 
v m (X - ray) 6 • 10 I6 Hz and v m (Optic) =s 2 ■ 10 14 Hz (where v m 
is the frequency of maximum in synchrotron emission spectrum) 
enter the same domains as the measured spectra. 

The original waves, excited during the cyclotron resonance, 
come in the radio domain, but the radio emission is not observed 
from RXJ1856. One of the possible explanations why the ra- 
dio emission is not detected from this object is that it traverses 
a large distance in the pulsar magnetosphere (since the model of 
a nearly aligned rotator is used). So there is a high probability 
for the excited waves to com e in the cyclotron damping range : 
to -kipVip - k x u x - u) B /y r = dKhechinashvili & Melikidze 111 9971) . 
In this case, the radio emission will not reach an observer. Never- 
theless, the detection of radio emission from RXJ1856 would be a 
strong argument in favour of the model. 

The effectiveness of the cyclotron mechanism has been esti- 
mated and it appears to be quite efficient. For effective generation 
of waves it is essential that the time during which the particles give 
energy to waves should be more than 1 /T (where F is the growth 
rate of instability). Generated radio waves propagate practically in 
straight lines, whereas the dipolar magnetic field lines deviate from 
its initial direction, and the angle 9 = k\\jk^ grows. 9 is the angle 
between the wave line and the line of dipole magnetic field. On the 
other hand, the resonance condition (1) imposes limitations on 9 
i.e. particles can resonate with the waves propagating in a limited 
range of angles. The estimations show that the following condition 
p > 3 ■ 10 9 cm should be fulfilled. As the instability develops at 



3 



the distances r ~ 10 9 cm, it follows that the excited waves lie in the 
resonant region long enough for particles to acquire pitch angles 
and generate the observed radiation. 

The total energy available for the conversion into pulsar emis- 
sion has been estimated, which approximately equals to E 
n\, a KR 2 pc j\,mc i . Here, n& is the Goldreich- Julian density at the pul- 
sar surface and R pc is the radius of the polar cap. The estimations 
show E a 5 ■ 10 32 erg/s, which is enough to explain the observed 
X-ray luminosity of RXJ1856. 

The recently discovered 7s pulsations of the X-ray emission of 
RXJ1856 has been explained in the framework of drift wave driven 
model. The main feature of this model is that the spin period of pul- 
sar might differ from the observable period (for RXJ1856 the real 
spin period is estimated to be ~ Is), as a consequence of the exis- 
tence of very low frequency drift waves in the region of generation 
of the pulsar emission. These particular waves are not detected but 
only result in a periodical change of curvature of the magnetic field 
lines, which in turn cause the change of observed radiation with a 
period of the drift wave. 

In the present paper the emission polarization of RX J 1856.5- 
3754 is investigated, as we suppose that the measurement of the 
polarization properties of RXJ1856 will make clear its real emis- 
sion nature. It was shown that both the X-ray and optical emissions 
are linearly polarized. The degrees of polarization for the X-ray 
emission approximately equals IT^ w 8% in the energy interval of 
(0.15 - 0.26)keV and II; * (84 - 79)% in the energy interval of 
(0.26 - 0.9)keV. The calculation of degree of polarization for op- 
tical emission yields 11^ « 86%. I f the emissi on of this source has 
a thermal nature then according to lHo I d2007 ) (this model, among 
other thermal emission models, gives the best match of the entire 
spectrum) the polarization of X-ray emission should essentially be 
equal to 100%. But if the emission of this source is generated by 
the synchrotron mechanism, it is expected that both the X-ray and 
the optical emissions will be linearly polarized with the frequency 
dependent polarization degrees, giving the values from a few per- 
cents up to 84%. 
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